Electron correlation effects and ferromagnetism in Iron 
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Electron correlation effects in Fe are analyzed using a first principles LC AO-scheme. In our ap- 
proach, we first use a local orbital DFT-LDA solution to introduce a Hubbard Hamiltonian without 
fitting parameters. In a second step, we introduce a many-body solution to this Hamiltonian using a 
DMF approximation. Our analysis shows that magnetism in Fe is an effect associated with the first 
atomic Hund's rule. Moreover, we also find important correlation effects in the Fe-spin polarized 
DOS. The photoemision spectra is explained using a value of C/ off as large as 4 eV, provided the 
satellite peaks appearing around 3-5 eV below the Fermi energy are interpreted appropriately. 
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The electronic properties of ferromagnetic metals are 
still a subject of controversy Although DFT- 

LDA calculations yield the correct magnetization for the 
itinerant-electron ferromagnets Fe, Co and Ni, the ori- 
gin of ferromagnetism in these metals and the role of 
electron correlations are not completely well-understood 
(e.g. see Q); in particular, the relative importance of the 
local Coulomb interaction for d orbitals, U eS , versus in- 
traatomic exchange (first Hund's rule) is not completely 
established. Even from an experimental point of view, 
there is a lack of agreement on the existence or not of 
a satellite peak in the photoemission spectrum of iron 
around 5 eV below the Fermi energy [9 . 

In the conventional view of itinerant ferromagnetism 
|^| , spin polarization is determined by the Stoner param- 
eter, /, that defines the energy of the atomic d-orbitals 
as (Ed — Irid), rid being the occupation number of the 
orbital under consideration. In the case of Fe, DFT- 
LDA calculations yield a value I — 3.9 eV M, and a 
surprisingly large value of C/ cff ~ 4-6 eV |7],|8| for the 
effective Coulomb interaction between the d-electrons. 
Since the atomic-like properties of the d-states are of 
crucial importance for the magnetic properties of these 
materials, LCAO methods provide the appropriate con- 
ceptual framework to understand those properties and 
to analyze the role of electron correlations. In LC AO- 
theories of ferromagnetism, including Hubbard Hamilto- 
nians, I is written as (U cS + 4J X ) [0], where J x defines 
the screened intrasite exchange interaction between the 
atomic d-electrons having the same spin; it is commonly 
accepted that J x practically coincides with its atomic 
value 0. In the case of Fe, J x = 0.83 eV and, therefore, 
we should take U cS ~ 0.6 eV to recover the value I = 3.9 
eV that corresponds to the correct magnetization. This 
result suggests the presence of dramatic electron corre- 
lation effects in Fe, which would be responsible for the 
renormalization of U from ~ 5 eV to U ~ 0.6 eV. 
On the other hand, the value of U cS inferred from the 
photoemission spectra [[H)]-[l2| by identifying photoemis- 
sion peaks with quasi-particle peaks yields U cS « 2 eV 
This result seems to indicate that electron correla- 



tion effects for Fe are not strong, in contradiction with 
the previous analysis. 

The purpose of this paper is to show that these ap- 
parent contradictions disappear once electron correlation 
effects fl3|| are properly analyzed using a first-principles 
LCAO-scheme. In our approach, reminiscent of LDA+U 
, we first formulate a Local Density solution for a gen- 
eralized Hubbard Hamiltonian. This LD-solution pro- 
vides the link between the generalized Hubbard Hamil- 
tonian and local orbital DFT-LDA methods, allowing us 
to obtain that Hamiltonian from first principles, with- 
out having to introduce fitting parameters. In a second 
step, we introduce a many-body solution for the Hub- 
bard Hamiltonian using a Dynamical Mean Field (DMF) 
approximation Jl5| : in this way we analyze the spin- 
polarized electron density of states for Fe and compare it 
with the experimental evidence JTc|— . From our anal- 
ysis, we obtain two different results. First, using our 
LD-solution for the Hubbard Hamiltonian, we show that 
electron correlation effects screen strongly the effective 
Coulomb interaction contributing to the Stoner parame- 
ter: in this sense, U is not larger than 0.6-0.7 eV. We 
find, however, that the effective interaction appearing in 
the Hubbard Hamiltonian is around 4 eV, in reasonable 
agreement with other first-principles calculations jjJD; 
using this value and the many-body techniques men- 
tioned above, we also find that the spin-polarized DOS 
for Fe is in good agreement with the photoemission data, 
provided we interpret appropriately the satellite peaks 
appearing in the spectrum around 3-5 eV below the Fermi 
energy [|. 

Our starting point is the generalized Hubbard Hamil- 
tonian: 



H = H OE + - Uih iaa h ll3al 
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where H defines a one-electron contribution, and Ui 
and Jiajp the intrasite and intersite coulomb interac- 
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tions between different orbitals <f>i a and (f)jp (for the sake 
of simplicity, Ui is an average of the different interactions 
inside the i-site); we also introduce the intrasite exchange 
coulomb interaction, Jf , associated with the first atomic 
Hund's rule. Eqn. |l| has been written in an orthogo- 
nal local basis, (f>i a , defined by the Lowdin's transforma- 
tion 4>i a — 'J2j l3 (S~ 1 / 2 )iajf3ipji3, ^jff being the local basis 
used in the DFT-LDA calculation from which we obtain 
Hamiltonian [l], as explained below (Si a jp is the overlap 
between orbitals ipia and ipjp)- 

The LD-solution of Hamiltonian |l| is obtained by intro- 
ducing the kinetic and many-body energies of the system 
as a function of the orbital occupancies, ni aa ||l6|| . This 
implies that the total energy is a function of those num- 
bers, ni aa , that play the role of the electron density, p(r), 
in the conventional DFT-approach. Then, we can write 
the following eqn.: 

E[{n iaa }] = T[{n iaa }] + E H [{n iaa }] + E xc [{n iaa }}, 

(2) 

where T = (^ \H OE \^ ), * being the ground state 
of the total LD Hamiltonian; E H is the hartree energy 
and E xc is the exchange-correlation energy associated 
with Hamiltonian [|. On the other hand jl6) : 



xc _ dE xc [{n laa }} _ v 
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an eqn. that yields the exchange energy as the sum of an 
intrasite contribution and of the intersite interaction be- 
tween the electron charge, rii ar7 , and its hole, (1 — rii aa ). 
In this eqn., Ji is practically the coulomb interaction be- 
tween charges located in n.n. atoms. Because of the crys- 
tal symmetry , we assume that no exchange hole appears 
in the same atom where the electron is located. On the 
other hand, we have also shown jl6| that the correlation 
energy is given by: 

E c [{n iaa }] = --^/i(£7j - Ji)n iaa (l - niaa) (4) 

IOLO 

where (Ui — Ji) is an effective intrasite coulomb interac- 
tion between i-site orbitals, and fi a quantity ranging be- 
tween and 1 depending on the importance of the intra- 
site correlation effects (fi is 1 for large values of (Ui — Jij). 
Eqns. |^, |^ and |I] allow us to substitute Hamiltonian [j] 
for an effective Hamiltonian where, instead of the many- 
body terms, we introduce the local potentials (V^ a and 
given by: 

H dE H [{n taa }} _ 



dn 

E 
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where fi has been assumed to be constant. 

This is the main result of our LD-analysis and shows 
how to reduce the generalized Hubbard Hamiltonian, 
eqn. pi to an effective one-electron Hamiltonian, tak- 
ing into account all the many-body contributions. Con- 
versely, we can use this equivalence to go from a LD- 
solution to a generalized Hubbard Hamiltonian: Assume 
we solve the conventional DFT-LDA eqns for a given 
crystal (say, paramagnetic Fe) using a local orbital basis 
(as done in Fireball |l7|] or Siesta [pj| codes); then, we 
can substract from the one-electron levels associated with 
the orthogonalizcd orbitals, (f>i a , the potentials given by 
eqns. H and ^|. This difference defines Hamiltonian H OE 
in eqn. |l|, and allows us to introduce the Hubbard Hamil- 
tonian by means of the interactions U, J and J x . Notice 
that in this approach we have to calculate these inter- 
actions using the orthogonalizcd orbitals, (f>i a . In our 
actual calculations, we have employed the Fireball code 
for paramagnetic Fe and used the corresponding local or- 
bital basis. 

Ferromagnetic Fe has been analyzed in our LD- 
approach by looking for a magnetic solution where some 
charge is transferred between spins up and down. This 
implies that selfconsistent potentials, Vi a \ and Vi a ±, 
should appear for different spins, in such a way that: 

Vi a ^ = Ui(Ni — n ia ^) + ^ Jia,j(3nj(3a' 
30tr' 

- J ii N i-\ ~ n *al) - Ji(^ ~ niaf) 



-fi(Ui - Jj)(- - n la] ) 



(7) 



UiUipa' + ^ Jj 



(5) 



where Ni is the total charge in the d-orbitals for the 
i-site, while represents the total spin-up charge. No- 
tice how the terms contributing to Vi a \ correspond to 
the Hartree, the intraatomic exchange, the extraatomic 
exchange and the correlation contributions, respectively. 
Due to the magnetic polarization, we find changes in the 
many-body potential w.r.t. the paramagnetic solution. 
This yields: 

&V ia] = -(1 - fi)(Ui - Ji)5n ia i - J?(5Nii - 5n ia] ) (8) 

where the total charge, Ni, has been assumed to be con- 
stant and independent from the atomic magnetization. 
Eqn. ^| defines how the ia "f -level depends on the atomic 
polarization, Sn^. This quantity should be obtained 
selfconsistently by means of a band-structure calcula- 
tion whereby the atomic charges, 5fii a -\, are a function 
of SVi a i ■ These two conditions yield 8ni a -\ and the crys- 
tal magnetization. Eqn. || allows us, however, to calcu- 
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late directly the Stoner parameter, /, which we define as 
\5Vi a -[ / 8rii a ^\. Eqn. || yields the following result: 



/=(l-/ i )(C/ i -J i )+4Jf 



(9) 



assuming SN^ — 56rii a -[, as corresponds to d orbitals. In 
our calculations for Fe, we find Ut = 14.7 eV (taking into 
account atomic relaxation), p = 6 eV and Jf= 0.83 eV. 

A word of caution should be introduced here, because 
in our discussion we have neglected an effect that leads 
to a further reduction in the effective interaction between 
orbitals. This is associated with the sp-band screen- 
ing that has been shown by other authors [7] to reduce 
(Ui — Ji) = U° s to values close to 5 eV. In our cal- 
culations, performed introducing a Lindhard dielectric 
function, we have found that J7f ff is reduced to 4.0 eV. 
If we introduce this value in eqn. | and take j l =0.83 
(the value that corresponds to this reduced interaction), 
we find / = 4.0 eV, in good agreement with DFT-LDA 
calculations. We have also analyzed how / depends on 
P° ff by calculating ji for different intrasite coulomb in- 
teractions. Our results show that, in the 2-5 eV range, 
(1 — fi)Ui S is almost insensitive to the values of Uf s . 
These results show that the Stoner parameter is mainly 
controlled by Jf; in other words, ferromagnetism in Fc 
is an effect mainly associated with the intraatomic first 
Hund's rule. 

Next, we calculate many-body effects introducing a lo- 
cal selfenergy, £i QO -(c<j), within the DMF-approximation. 
This is a reasonable approximation considering that cor- 
relation effects in Fe are associated with the intrasite 
coulomb interaction between d orbitals. As discussed in 
ref 1 16 1, Tii aa {ijj) is calculated by an appropriate inter- 
polation between two limits: (1) first , we calculate the 
atomic limit, assuming Uf s much larger than the metal 
bandwidth; (2) second, we obtain the second order selfen- 
ergy, S-^^w), using as the expansion parameter U° s ; (3) 
finally, we calculate the selfenergy interpolating between 
these two limits. 



We should stress that, in this solution, £ 



(2) 



is cal- 



culated using the local density of states defined by the 
LD-solution discussed above. Finally, we replace Vf aa , in 
our effective LD-Hamiltonian, by that selfenergy; then, 
we use a conventional Green-function formalism to calcu- 
late the local density of states. At this point, we should 
comment that consistency between the LD- and the self- 
energy formulations imposes the following Luttinger sum 
rule: Vf a(J — Y, iaa (Ep). The factor /, in eqn ||] has been 
determined from this eqn., this procedure satisfying the 
Luttinger condition automatically. The price we have 
to pay is the introduction of a selfconsistent loop in the 
calculation. 

Fig. 1 shows our calculated LD- and many-body DOS 
for ferromagnetic Fe. The comparison between these two 
density of states shows that correlation effects are im- 
portant for Fe: first, we notice that the energy difference 



between the two maxima appearing in the spin-up and 
spin-down DOS for the LD-solution is reduced by almost 
a factor of two in the many-body case. This is a typical 
band narrowing effect appearing around Ep and associ- 
ated with a highly correlated electron gas. On the other 
hand, we also find that the DOS-structure at energies far 
away from Ep is strongly modified by the many-body 
solution: in particular, a new satellite structure appears 
around 5 eV below Ep. We analyze more in detail these 
many-body effects by considering the DOS for a partic- 
ular fc-vector: we have chosen the P-point, a case for 
which there are high quality photoemission data taken 
along the (111) direction 0Q. 
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FIG. 1. Total Fe spin DOS for the many-body solution with 
U cB = 4 eV (continuous line) and for the DFT-LDA solution 
(dashed line). Notice the shift of the peaks towards Ef when 
correlation effects are introduced. 




++++- 
-4 -2 
co(eV) 

FIG. 2. Spin DOS, split into the t2 g and e g components, at 
the P-point for Fe. The position of the DFT-LDA eigenvalues 
is indicated by arrows. Notice that the t2 ff (up) state (black 
squares) has lost its quasiparticle character. 



Fig. 2 shows the electron DOS for this fc-vector split 
^-components, and the energy levels 



into the ti g and the e 



corresponding to the DFT-LDA-solution for this /c-point. 
In this figure we also find the effects already discussed 
for fig. 1: the e g (up) and t2g(down) levels located 0.6 
eV and 2.1 eV below Ep are shifted to around 0.4 and 
1.4 eV, respectively. In both cases, we also find some 
satellites at higher binding energies due to many-body 
effects. The other £2g(up)4evel located 3.1 eV below Ep 
is completely smeared out by the selfenergy, giving rise 
to two features: one is almost coinciding in energy with 
the t2g(down)-peak, the other one is a very broad peak 
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located around 4.5 eV below Ep. For energies above 
Ep, we find the e g (down)-level, located 1.4 eV above Ep, 
shifted by many-body effects to 0.9 eV. In fig. 3 we com- 
pare these results with photoemission data |ll],[f^], by 
considering the appropriate weight that each state has 
in the photoemitted electrons: this is done assuming the 
final state to be a plane wave; in this way we find that, 
for the P-point, the i2g-levels are reduced in intensity by 
a factor of 3 w.r.t. the e s -levels. Comparing this DOS 
with the photoemitted spectra, we clearly see that the 
two peaks below Ep (at 0.4 and 1.4 eV) are related to 
the e ff (up) and (down) states calculated in LD, while 
the peak above Ep (at 0.9 eV) corresponds to a e g (down) 
level: these results are in reasonable agreement with the 
experimental data of refs. More importantly, we 

find that the very broad peak located around 3-4 eV be- 
low Ep cannot be related directly to the i2g(up)-states 
found in LD around 3.1 eV below Ep [jll|. On the con- 
trary, our results clearly show that this peak is a satellite 
structure created by many-body effects and appearing as 
the result of combining the tail intensities of e g and ti g 
states (see fig 2, and remember the factor 1/3 we have to 
introduce in the weight of the <2g -states). 
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FIG. 3. Quasiparticle spectra (considering the appropriate 
weight for each state, see text) at the P-point for Fe. Black 
arrows indicate the eigenstates of the DFT-LDA Hamiltonian. 
Experiments p""Ij|il^ show three well-defined peaks close to Ef 
and a very broad feature around 3 eV below Ep traditionally 
interpreted as the £23 (up) state (Gray arrows). 

This analysis clarifies several contradictory points 
about the magnetism of iron. Our results show that we 
can find a reasonable agreement between theory and pho- 
toemision data for U cS as large as 4 eV. The reason why 
a value of U cS ~ 2 eV has been used in the interpretation 
of these data is the tendency to identify the broad satel- 
lite peak, located around 3 eV below Ep, at the P-point, 
with the quasiparticle level that in DFT-LDA appears 
around 3.1 eV below Ep. Our results show, however, 
that this quasiparticle level has lost its identity due to 
correlation effects and that it has been modified into a 
smeared DOS: this new DOS tends to create a satellite 



structure that should be reinterpreted as due to many- 
body effects and not as a quasiparticle level reminiscent 
of the DFT-LDA level. This same effect is also responsi- 
ble of the satellite peak we find in the total DOS (see fig 
1) around 5 eV below Ep. 

In conclusion, we have studied the electronic proper- 
ties of ferromagnetic Fe using a first-principles LCAO- 
scheme to analyze in detail the role of electron corre- 
lations. We find that the correlation potential strongly 
screens the magnetic effects commonly associated with a 
local Hubbard interaction: magnetism in Fe is an effect 
associated with the first atomic Hund's rule. Moreover, 
we also find important correlation effects in the spin- 
polarized DOS. In particular, our analysis shows that 
the <2g(up) level located 3.1 eV below Ep have lost their 
quasiparticle identity, due to many-body effects, and tend 
to create a satellite structure that has been observed ex- 
perimentally. 
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CYT under contract No. PB-97-0028. F. F. thanks N. 
H. March for an illuminating discussion. 
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